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Variational wave functions are very useful for describing the panoply of ground states found in 
interacting many-electron systems. Some particular trial states are "adiabatically" linked to a refer- 
ence state, from which they borrow the essential properties. A prominent example is the Gutzwiller 
ansatz, where one starts with the filled Fermi sea. A simple soluble example, the anisotropic XY 
chain, illustrates the adiabatic continuity of this class of wave functions. To describe symmetry 
breaking, one has to modify the reference state accordingly. Alternatively, a quantum phase transi- 
tion can be described by a pair of variational wave functions, starting at weak and strong coupling, 
respectively. 



I. INTRODUCTION 



o 

(N 
o 

, According to the adiabatic theorem of quantum mechanics a system governed by a Hamiltonian with time-dependent 
00 ' parameters evolves from an initial eigenstate to an eigenstate at later times if the parameters vary infinitesimally slowly 
\ with time 0] ■ An adiabatic evolution from an initial non-interacting eigenstate to an interacting state by slowly 

■ switching on the interactions is also a key element in Landau's theory of the Fermi liquid [s', Anderson generalized 
1— H '■ this idea to the concept of adiabatic continuity, which basically means the continuous mapping of a simple reference 

^ . state - the non-interacting Fermi sea in Landau's theory - to a more complex state which however shares the essential 
' properties with the parent state 5] . This can only work for systems which do not experience any symmetry breaking 
c/3 , on the way. It is the aim of this paper to show that a similar continuity characterizes a whole class of variational 
wave functions. When these are applied without precaution to systems that undergo some phase transition one may 
easily be led to wrong conclusions. 

Section 2 recalls the well-known example of the Gutzwiller wave function, which represents a metallic state, even 
^ ' for a half-filled band. It can therefore be viewed as a microscopic realization of Landau's phenomenological theory 
(— ( , of the Fermi liquid. Section 3 presents the exactly solvable anisotropic XY chain which exhibits a quantum phase 
Q ' transition for an infinitesimal anisotropy. A Gutzwiller-type variational ansatz fails in signaling this transition. Section 
O . 4 deals with various kinds of quantum phase transitions, both with and without symmetry breaking. "Dual" pairs of 
variational wave functions give an appealing picture of metal-insulator transitions. In Section 5 a different class of 
I \ wave functions is briefly mentioned, which is able to promote long-range order without forcing symmetry breaking in 
J> ' the reference state. 
00 

■ II. ADIABATIC CONTINUITY OF THE GUTZWILLER ANSATZ 

When Martin Gutzwiller introduced his famous ansatz for the ground state of a correlated electron system he 
wanted "to present a new approach to the problem of ferromagnetism in a metal" Q. As a particular example he 
considered the Hubbard Hamiltonian 

H = —t ^ (cJj^Cmo- -I- cl^^Cna) + ^ Cj^^C„-|-c]^j^C„4, , (1) 
{n,m),(T n 



where cj^g. (ca) creates (annihilates) an electron with spin a at site n . The issue was not to find a good trial state for 
ferromagnetism, because the lowest-energy fully polarized ferromagnetic state for N electrons is simply constructed 
by occupying the N lowest levels with only up (or down) electrons, but rather to show that no other state has lower 
energy. This problem has been intensively studied ever since Gutzwiller's work and is still not completely understood 
0. Only for special lattices and/or band structures it can be rigorously shown that the ground state of the Hubbard 
model is ferromagnetic Q. 

Gutzwiller proposed the ansatz 



=exp 



l*o) (2) 
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for the correlated non-magnetic state, where |^'o) is the filled Fermi sea. The variational parameter g reduces double 
occupancy and thus the interaction energy. As to the kinetic energy, a simple approximate expression can be obtained 
by neglecting the dependence of the hopping processes on the detailed spin configurations [1, Q . This "Gutzwiller 
approximation" allows the calculation of physical quantities such as the momentum distribution and, by including a 
Zeeman term, the magnetic susceptibility. For a generic band filling, i.e., away from half filling, correlation effects 
essentially increase the effective mass but do not destroy the metallicity. An interesting link to Landau's Fermi-liquid 
theory has been established by VoUhardt fl^, who was able to establish a direct connection between microscopic 
quantities and Landau parameters. The smooth connection between the correlated state \^{g)) and the uncorrelated 
state Y^q) in Eq. ^ has thus its counterpart in the adiabatic hypothesis underlying Landau's Fermi-liquid theory. 

The special density of one electron per site (half-filled band) is worth being discussed in some detail, because in this 
case the Gutzwiller approximation predicts a metal- insulator transition at a critical value Uc = 8|eo|, where eo is the 
kinetic energy per site of the uncorrelated system [llj . Signatures of the transition are both the vanishing of double 
occupancy and the disappearance of the Fermi step in the momentum distribution function. However, more accurate 
treatments of the Gutzwiller ansatz obtain this transition only for ?7 — > oo and thus show that the variational state ^ 
remains metallic for all finite values of U . This is clearly seen in the case of the one-dimensional Hubbard model for 
which the Gutzwiller ansatz can be handled exactly 1^, 13 1, but there exists also a simple argument why this remains 
true for any (finite) dimension. The argument is based on the Drude weight which is related to the sensitivity of the 
ground state with respect to changes in the boundary conditions p^ . The Drude weight is finite for the Gutzwiller 
ground state for any finite value of U, even at half filling [15|, (l^ . 



III. QUANTUM PHASE TRANSITION FOR INFINITESIMAL COUPLING: THE ANISOTROPIC XY 

CHAIN 

A simple example illustrates the smooth link between a variational wave function of the Gutzwiller type and its 
reference state, the "bare" ground state. We consider the anisotropic XY chain, represented by the Hamiltonian 

H = Ho+-fH', (3) 

where 

Ho = Y.{sk^^Snll+Si'^S^"ll) ' 

n 

H' = y: - siy^sii,) . (4) 

n 

The parameter 7, < 7 < 1, interpolates between the isotropic XY chain (7 = 0) and the one-dimensional Ising 
model (7 = 1). This Hamiltonian has been diagonalized exactly by Lieb, Scliultz and Mattis The main steps 

of the solution are reproduced in the appendix. It turns out that this model exhibits a quantum phase transition at 
7 = 0, from a ground state without long-r ang e order at 7 = to a state with long-range antiferromagnetic order of the 
x-components of the spins for any 7 > 17]. The transition reveals itself also in the singularity of the ground-state 
energy (per site, 7 <C 1) 
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(5) 



Another interesting quantity is the overlap (\I'(7')|4'(7)) between ground states of two Hamiltonians with different 
anisotropics [l^. We consider in particular the overlap at criticality, i.e. for 7' = and 7^1. It is convenient to 
introduce the fidelity susceptibility per site [Toj 



XF= limT7^1og(M'(0)|vI'(7)). (6) 



-2 

7->-6 Lj 

As shown in the appendix, we find for the XY chain for 7 1 



which shows that xf diverges, as expected for a quantum-critical point [20l |. 
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We turn now to a variational ansatz a la Gutzwiller. The transformation to fermion operators c„, cj^ (appendix) 
leads to the expressions 

n 

H' = ^^(44+1 +c„+ic„), (8) 

n 

or, after Fourier transformation, 
Hq = ^ cos k c\ck , 

k 

H' = sin fc (cJ.cLfc - c^fcCfc) . (9) 

fc>0 

The Gutzwiller ansatz 

|vI/(5)) = e-f^>o> , (10) 

where |5'o) is the ground state of Hq, i.e., all levels with > tt/2 are occupied and all others are empty, can be 
written 

1^(5))= n (coshi9fc -isinhi9fe4cLfc) (coshi?fe + isinhi9fcC_fcCfc) l^-o) , (11) 

0<fc<7r/2 k>w/2 

where -dk — gsink. The expectation value of the Hamiltonian is readily calculated, 

{^(g)\H\^(g)) V- f Icosfcl . , • ,n1 

- E .Jr.„J.,, +7smfctanh(2.gsmfc) . (12) 



(*(5)l*(3)> ^ lcosh(25sinfc) 

The minimization of this expression for g <C 1 gives g = (37r/8)7 and an energy per site 

This is an analytic function of 7, in contrast to the exact result, Eq. ([S]). A striking difference is also obtained for the 
overlap. The (normalized) trial state of Eq. PT|) leads to 

immg)) -TT cosh(gsmfc) 1 
v/(*(5)l*(3)> v/cosh(2gsinfc) 2^ ^ 

Thus in the limit 7—7-0 the fidelity susceptibility tends to a constant for the Gutzwiller ansatz. 



XF^l^Y^j ~ 0.35, (15) 
while it diverges ^7^^ for the exact ground state. 

IV. VARIATIONAL WAVE FUNCTIONS AND QUANTUM PHASE TRANSITIONS 

The discrepancy found for the anisotropic XY chain between the exact and variational ground states exists also in 
the case of the one-dimensional Hubbard model. At half filling the exact ground state energy [2l| is not analytic for 



[/ — )• |22j, J while the Gutzwiller ansatz yields an analytic behavior [2J| . It is natural to attribute this difference 
to the fact that a metal-insulator transition occurs in the exact ground state at J7 = 0, but not so in the Gutzwiller 
ansatz, as explained above. We can remedy this problem by modifying the reference state |vl/o), which so far was 
considered to be metallic. We could introduce an alternating spin density, as predicted by the Unrestricted Hartree- 
Fock approximation. A more interesting possibility is bond alternation in the sense of an alternating bond order. This 
instability does not occur in a simple mean- field treatment, which would require a finite electron-phonon coupling 
(Peierls instability). However, if the reference state \'i>o) is taken as a mean- field state with alternating bond order. 
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the (Gutzwiller) variational energy is lowered, even in the absence of electron-phonon coupling. The condensation 
energy due to this symmetry breaking is given by the non- analytic term [251 ] 

£co„<i~-ie-3-«^(4*/^)\ (16) 

Long-range order is of course not expected to survive if quantum fluctuations are fully taken into account, rather one 
expects the correlation functions, both for spin and for the bond order, to decay algebraically with distance. 

In two dimensions fluctuations are less severe and therefore long-range antiferromagnetic order, as found consistently 
for the two-dimensional Hubbard model at half filling by using simple mean-field theory or more elaborate methods, is 
commonly believed to resist quantum fluctuations. A more subtle question is the issue of superconductivity because, 
similarly to bond alternation in one dimension, a superconducting instability does not occur within simple mean-field 
theory. For certain electron densities the Gutzwiller ansatz, linked to a reference state with d-wave pairing, is found 
to have lower energy than the state without symmetry breaking (or with an antiferromagnetic reference state) [26| . 
This remains true for more sophisticated variational ground states [27l - [29| . Whether the exact ground state exhibits 
superconductivity remains an open issue (soj . but the variational results indicate at least strong superconducting 
correlations. 

We turn now to the Mott phenomenon, which a priori does not involve any symmetry breaking, but rather is a 
topological (quantum) phase transition as a function of coupling strength, from a metallic to an insulating state. The 
transition typically occurs in a region where the interaction strength is of the order of the band width. In contrast 
to the instabilities discussed above, the Mott metal-insulator transition cannot be described in the framework of the 
Gutzwiller ansatz alone, which is metallic in the absence of symmetry breaking, but rather calls for a complementary 
variational state for the insulating side of the transition. In the case of the Hubbard model the following ansatz has 
been proposed for the ground state of the insulating phase [sTj , 



\^{h)) = exp 



(n,m) ,cr 



l*oo), (17) 



where |5'oo) is the exact ground state for U — > co. This state, the dual partner of the Gutzwiller ansatz [321, 
also adiabatically connected to its reference state |^oo)- The Drude weight vanishes at half filling [16] and therefore 
\^(K)) does represent an insulating state. Unfortunately, it is very hard to handle l^'oo), which is the ground state of 
the Heisenberg antiferromagnet at half filling and of the t-J Hamiltonian close to half filling, except for some special 
cases, such as the 1/r Hubbard chain "s^] or for infinite dimensions [l6j . In these cases a variational procedure using 
two trial states, one linked to the ground state at [/ = 0, the other to the ground state for U — oo, is rather successful 
in predicting the location of the Mott transition [l^. On the other hand, the nature of the transition is inevitably 
found to be of first order, also in cases where it is known to be continuous [32| . 

A related phenomenon is the crystallization of electrons, proposed by Wigner for the dilute three-dimensional 
electron gas [3j] and later by Hubbard for narrow-band solids '34!| . In this case the metal- insulator transition leading 
to the Wigner crystal involves a reduction of translational symmetry. For the case of electrons on a lattice the trial state 
of Eq. ((T7)) can again be used, this time for describing the crystalline phase, starting from the classical configuration 
which minimizes the Coulomb energy [ss'l . This leads to a quantum crystal where electrons are preferentially located 
at the sites of the reference configuration, but occasionally visit neighboring sites. With increasing relative strength of 
the kinetic energy these excursions become more and more frequent until the crystal melts. It would be interesting to 
test this two-sided variational procedure for the anisotropic XY chain by studying the dual partner of the Gutzwiller 
ansatz, starting from the Neel state, the ground state of the antiferromagnetic Ising model. 



V. DISCUSSION 



The variational ground states presented above are all smoothly linked to some reference state. The case of a tran- 
sition occurring at a large value of the coupling constant, such as the Mott transition and the Wigner crystallization, 
is rather well described by a pair of wave functions that are linked to the ground states for vanishing and infinite 
coupling strengths, respectively. In this way the critical coupling strength may be determined rather accurately, 
but the nature of the transition is not described reliably. Therefore more refined methods are required for treating 
the region of the transition, which may be of second order, in contrast to the prediction in terms of two dual wave 
functions, or it may even be split into two critical points with an intermediate phase in-between, as debated in the 
case of Wigner crystallization |36l . 137^ . 

A single variational state of the type discussed above does not change the qualitative nature of a reference state. 
For instance, these states do not lead to long-range order not included already in the reference state |^'o), nor do 
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they destroy long-range order if it is present in \^o)- To achieve such a radical change due to a continuous change of 
some variational parameter, a different class of variational states has to be used, such as the resonance-valence-bond 
state with a variable bond-length distribution [s^- This state is able to describe both short-range and long-range 
antiferromagnetic order on a square lattice. It is also possible to describe the Mott transition using a single variational 
wave function, by adding a Jastrow factor correlating distant electron densities [39|] • 
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Appendix A: Anisotropic XY model in one dimension 

We consider the anisotropic XY chain defined by the Hamiltonian 

L 

^ = E [(1 + + (1 - 7)5(^^5^1] , (Al) 

n=l 

where si"^ are spin 1/2 operators and < 7 < 1. We assume periodic boundary conditions, Si'L — S"}" , 

a = x,y, and choose L — AM + 2 (where M is a positive integer). The Jordan- Wigner transformation 

^(-)-h^5^«) = cte"S":^^^ (A2) 

where cjj, c„ are fermionic creation and annihilation operators, yields a simple quadratic expression for the Hamilto- 
nian, 

1 ^ 

^ = 2 E(^«^"+i + ^44+1 + h.c.) . (A3) 

n=l 

Introducing the Fourier transform 

c„ = ^ y e^'^^Cfc, k = ^i/, -2M <i^<2M + l, (A4) 
we rewrite the Hamiltonian 

= y cos fc c[cfc + 17 E ^ i^k^-k ~ c^kCk) (A5) 

k k>0 

and diagonalize it using the Bogoliubov transformation 

Cfc e*''/'*(cos ipk dk + sin i^k c^Lfc) , 
cik = e"*''/''(-sin(^fcdfc -|-cos(pfcdLfc), (A6) 
where 

tan(2(y5fe) — — 7tanfc. (A7) 
This yields the excitation spectrum 



Ek = ^\- (l-72)sin2A:. (A8) 
With the choice 

'^'^'^=[ 2E, ) ' ""^^ = '1, 2E, ) 
the Hamiltonian reads 

H = J2Ek(didk-l) . (AlO) 
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The ground state energy per site for L — >■ oo is given by 



(All) 



where E{1 — 7^) is the complete elliptic integral. The ground state [^'(7)), defined by dk\'^{'~f)) = , can be written 



(A12) 



fe>0 



where |0) is the vacuum of c-particles. The overlap between two ground states for different anisotropies 7,7' is given 
by 



(M'(7')|*(7)> = n '^^^(^^ - ^fc) • 



(A13) 



fc>0 



We are particularly interested in the overlap between the ground state of the isotropic XY model (7 — 0) and that of 
the anisotropic model (7 > 0). With the choice (jA9p we get 



(*(o)i*(7)) = n 



k>0 



Ek + I cosfc|^ 2 
2Ek 



i.e., in the thermodynamic limit, 



log(vI/(0)|vI/(7)) 



L 

2^ Jo 



dk log ■ 



Ek + cos k 



2Ek 



(A14) 



(A15) 



To obtain the asymptotic behavior for 7 — >■ we have to distinguish between the region of k values close to ^ and the 
other region in k space. In the latter we can readily expand the integrand in powers of 7, while in the former region 
we replace cos k hy ^ ~ k and then integrate that part exactly. Summing up the two contributions we find 



log(vl/(0)|*(7)> ^ ^ 



-77-1 7+777 



(A16) 
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